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ABSTRACT 

We present the results of a photometric survey of rotation rates in the Coma Berenices 
OO (Melotte 111) open cluster, using data obtained as part of the SuperWASP exoplane- 

tary transit-search programme. The goal of the Coma survey was to measure precise 

rotation periods for main-sequence F, G and K dwarfs in this intermediate-age 600 
QQ Myr) cluster, and to determine the extent to which magnetic braking has caused the 

stellar spin periods to converge. We find a tight, almost linear relationship between 

rotation period and J — K colour with a root-mean square scatter of only 2 percent. 

The relation is similar to that seen among F, G and K stars in the Hyades. Such strong 
^ convergence can only be explained if angular momentum is not at present being trans- 

• ^ ferred from a reservoir in the deep stellar interiors to the surface layers. We conclude 

that the coupling timescale for angular momentum transport from a rapidly-spinning 
J-j radiative core to the outer convective zone must be substantially shorter than the 



cluster age, and that from the age of Coma onward, stars rotate effectively as solid 
bodies. The existence of a tight relationship between stellar mass and rotation period 
at a given age supports the use of stellar rotation period as an age indicator in F, 
G and K stars of Hyades age and older. We demonstrate that individual stellar ages 
can be determined within the Coma population with an internal precision of order 
9 percent (RMS), using a standard magnetic braking law in which rotation period 
increases with the square root of stellar age. We find that a slight modification to the 
magnetic-braking power law, P cx t"-^^, yields rotational and asteroseismological ages 
in good agreement for the Sun and other stars of solar age for which p-mode studies 
and photometric rotation periods have been published. 

Key words: stars: rotation - stars: activity - open clusters and associations: indi- 
vidual: Melotte 111 - techniques: photometric - methods: data analysis 



1 INTRODUCTION 
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In their pioneering study of stellar rotation among main- 
sequence stars in the Hyades open cluster, [Radick et al.| 
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( 1987 1 found a surprisingly narrow correlation between the 
B — V colours and the rotation periods of F, G and K 
stars. This discovery has several far-reaching implications. 
First, it confirms the theoretical prediction that the spin 
rates of an ensemble of stars of the same age and mass, 
but different initial rotation periods rates, should converge. 
This comes about because the rate at which angular mo- 
mentum is lost through magnetic braking via a thermally- 
driven, magnetically-channelled wind, increases strongly as 
a function of stellar rotation rate (e.g. Mestel & Spruit 
1987). Radick et al showed that this convergence is essen- 
tially complete at the ~ 600 Myr age of the Hyades. Second, 
for convergence to occur at such a comparatively early age, 
the timescale for coupling the star's radiative interior to its 
outer convective zone must be short. Models that allow a 
star's outer convective zone to be spun down rapidly by 
wind losses at the 50-100 Myr ages of clusters such as the 
Pleiades, retaining a hidden reservoir of angular momentum 
in a rapidly-spinning radiative interior, do not show this con- 



M37, Hartman et al. ( 2008 1 identified a substantial clump 



vergence (Li fc Collier Cameron|[l993 Krishnamurthi et al 



1997 Allain||1998 1 



For the core and envelope to be decoupled at the age of 
the Pleiades but to be rotating synchronously by the age 
of the Sun, the coupling timescale must be no less than 
~ 10 Myr, and no more than 1 Gyr ( [Soderblom, Jones, 
\&c Fi scher"2001'). Rotational evolutionary tracks with cou- 
pling timescales of order 100 Myr do not converge, because 
the "buried" angular momentum re-emerges at ages of a few 
hundred Myr. The wide range of spin rates seen in the very 
youngest clusters then re-asserts itself, inducing a large scat- 
ter in the predicted rotation rates and delaying convergence 
well beyond the age of Coma and the Hyades. 



The Hyades study by Radick et al. ( 1987 1 remained un- 



parallelled for nearly two decades after its publication. The 
advent of wide-field synoptic surveys for exoplanetary tran- 
sits in open clusters using Im-class telescopes, and in nearby 
field stars using even smaller instruments, has led to a recent 
resurgence of interest in the rotation distributions of open 
clusters with ages of ~ 150 Myr and older. Several recent 
photometric studies show the convergence of spin rates in 
nearby open clusters of intermediate age, notably M35 (150 



Myr; Meibom, Mathieu, fc Stassun 20091, M34 (250 Myr 



Irwin et al.||2006|) an d M37 (550 Myr: [Hartman et al.||2008 



Messina et al. 20081. These studies have gained new sig- 



nificance because the rotation-age-colour relations in these 
clusters provide important calibration points for stellar ro- 
tational evolution. An empirical calibration of stellar spin 
rates as a function of age and colour was published recently 
by [Barnes (2007), who coined the term gyrochronology to de- 
scribe the inverse process of inferring a star's age from its ro- 
tation period and colour. If magnetic braking causes stellar 
spin rates to converge by the age of the Hyades (or perhaps 
sooner) and subsequently to spin down as some unique and 
measurable function of time, gyrochronology becomes very 
attractive as a clock for measuring the ages of individual 
field stars, stellar and planetary systems. 

While this convergence of spin rates seems to be essen- 
tially complete among F and G stars by the age of M35 
( Meibom, Mathieu, fc Stassun|2009 1, and in K stars by the 
age of the Hyades, it takes even longer for stars of still lower 
mass. Although many stars with B — V > 1.3 have converged 
to the main period-colour relation by the 550-Myr age of 



of stars redder than this limit having periods of 2.5 days or 
less. In their preliminary investigation of the still older Prae- 
sepe (650 Myr), |SchoIz"fc Eisloffel] (|2007[) find convergence 



to be complete down to the late K types, but rapid rotation 
persisting among low-mass stars of spectral type MO and 
later. 

Here we present the results of a wide-field photomet- 
ric monitoring survey of stars brighter than V — 14 in 
the region around the Coma Berenices open cluster, using 
the SuperWASP camera array on La Palma. This cluster, 
also known as Melotte 111, is centred at RA = 12h 23m, 
Dec. = 4-26° 00 (J2000.0). It is the second-closest open 
cluster to the Sun, lying at a distance of 89.9 ±2.1 pc 



van Leeuwen 1999), further only than the Hyades (46.3 
pc: [Ferryman et al.||1997| . The foreground reddening is 



E{B - V) ~ 0.006 ±0.013 ( |Nicolet[[l98l| ) . [Trumpler| ( |l938 ) 
used proper-motion, radial-velocity and colour-magnitude 
selection criteria to identify 37 probable members with pho- 
tographic magnitudes nipg < 10.5 within 3.5° of the cluster 



within the same region in a proper motion survey by |Ar-| 



tiukhina (1955), but only two of these were confirmed by 



Argue & Kenworthy ( 1969 ), who identified another two faint 



objects that were later confirmed as members. |Odenkirchen,| 
Soubiran, fc Colin (!l998) performed a thorough kinematic 
and photometric survey down to 1/ ~ 10.5 using the Hip- 
parcos and the combined Astrographic Catalogue and Tycho 



(ACT; Urban, Corbin, fc Wycoff 1998 ) reference catalogues. 
This study yielded around 50 probable kinematic members 
in a region of 1200 square degrees about the cluster cen- 
tre, effectively superseding the kinematic and photometric 
aspects of Trumpler's original survey. [Casewell, Jameson, fc| 



Dobbie ( 2006 1 compiled a list of previously known mem- 



bers (according to proper-motion, photometric and in some 
cases radial-velocity criteria) within 4 degrees of the cluster 
centre. They found a further 60 candidate cluster members 
using the USNO-Bl.O (United States Naval Observatory; 



Monet et al. 2003) and 2MASS (Two-Micron All-Sky Sur- 



vey) point-source catalogues to carry out proper motion and 
photometric surveys respectively, more than doubling the 
number of probable Melotte 111 cluster members to about 
110 according to photometric and kinematic criteria. Among 
these, about 40 are secure in the sense that they also satisfy 
radial-velocity membership criteria, at least at the level of 
precision achieved by Trumpler (1938). 



Previous photometric studies by Radick, Skiff, & Lock- 



wood ( 1990 ) of four early-G dwarfs and by Marilli, Cata 



lano, fc Frasca[ ( [1997[ ) of three late-G dwarfs in the Coma 



Berenices cluster suggest a Hyades-like convergence in spin 
rates. Like the Hyades, however. Coma Berenices is a diffi- 
cult target for high-precision photometric rotation studies, 
because of its age and and its wide angular extent. As in the 
Hyades, the stellar spin periods are expected to range from 
6 days among late F dwarfs to 14 days among mid-K dwarfs, 
with more rapid rotation persisting at later spectral types. 
Although stellar dynamos still produce strong magnetic ac- 
tivity signatures at such rotation rates, the starspot coverage 
is substantially lower than in the ultra-fast rotators found in 
younger open clusters. The comparatively low starspot cov- 
erage on stars in the Hyades and Coma yields amplitudes of 
rotational modulation of at most one or two percent. This 
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combination of long periods, low amplitudes and large angu- 
lar separation between cluster members presents a challenge 
to observers seeking precise measurements of their optical 
modulation periods. 

In sections 2 and 3 of this paper we describe the ob- 
servations and data reduction methodology. In section 4 we 
describe the generalised Lomb-Scargle period-search method 
we used to identify candidate rotational variable stars in the 
Coma Berenices region, while in section 5 we describe the 
additional photometric and astrometric criteria used to as- 
sess their cluster membership probabilities. In Sections 6 and 
7 we discuss the period-colour relation for Coma Berenices 
cluster members, and the implications for stellar spindown 
and gyrochronological age determination. 



2 OBSERVATIONS 

The SuperWASP camera array, located at the Observatorio 
del Roque de los Muchachos on La Palma, Canary Islands, 
carried out its inaugural season of observations during 2004. 
The camera array is mounted on a robotic equatorial mount. 
In 2004 it comprised five 200mm f/1.8 Canon lenses each 
with an Andor CCD array of 2048'^ 13.5/xm pixels, giving a 



Table 1. Journal of SuperWASP observations for the Coma 
Berenices rotation survey. 



field of view 7.8 degrees square for each camera (PoUacco et 
aL]|20 06'). Four individual camera fields covering the region 
around the Coma Cluster were observed on 51 nights of good 
photometric quality between 2004 May 2 and 2004 July 18; 
a further two fields at an adjacent pointing position were 
observed on 86 nights between 2004 May 02 and 2004 August 
12. 

The Coma Berenices region was re-observed throughout 
the first half of 2007. By this time a further three cameras 
had been installed on the mount. Again, six individual cam- 
era fields provided full coverage of the region around the 
cluster: two on 72 nights from 2006 December 12 to 2007 
May 15, and four from 2007 January 14 to May 30. 

The survey region was accessible for between 4 and 8 
hours each night. The average interval between visits to the 
field was 6 minutes. Each exposure was of 30s duration, and 
was taken without filters in 2004, and with an infrared block- 
ing fiher in 2007. 

The field centres and dates of observation are sum- 
marised in Table [l] 



3 DATA REDUCTION AND CALIBRATION 

The data were reduced using the SuperWASP data reduc- 



tion pipeline, whose operation is described in detail by (Pol 



lacco et al. 2006). Science frames are bias-subtracted and 
fiat-fielded. An automated field recognition algorithm iden- 
tifies the objects on the frame with their counterparts in the 
TYCHO-2 catalogue and establishes an astrometric solution 
with an RMS precision of 0.1 to 0.2 pixel. The pipeline then 
performs aperture photometry at the positions on each CCD 



image of all objects in the USNO-Bl.O catalogue (Monet et 



|al.|2003) with magnitudes in the R2 band brighter than 14.5, 
equivalent to V ~ 15 over the colour range < B — V < 1. 
The pipeline computes he formal variance of the flux of each 
star in the image from the photon counts in the stellar aper- 
ture and the surrounding sky annulus. This magnitude limit 
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is dictated by the detection threshold for the SuperWASP 
exoplanetary transit-search project; stars with R2 ~ 14.5 
yield light curves with RMS scatter between 0.02 and 0.03 
magnitude. 

The SuperWASP bandpass covers most of the optical 
spectrum, so the effective atmospheric extinction depends 
on stellar colour. The extinction correction to zero airmass 
for a star of known B — V colour index takes the form 

mo = m - k'X - k"{B - V)X 

where m is the raw instrumental magnitude at airmass X, 
mo is the instrumental magnitude above the atmosphere, 
and k' and k" are the primary and secondary extinction 
coefficients respectively. The zero-point z for each image is 
tied to a network of local secondary standards in each field. 
The resulting instrumental magnitudes are transformed to 
TYCHO-2 V magnitudes via the colour equation 



Vt = mo 



c{B - V) 



where the colour transformation coefficient for all cameras 
lies in the range c = —0.22 ± 0.03. The resulting fiuxes are 
stored in the SuperWASP Data Archive at the University of 
Leicester. 

We extracted the light-curves of all stars observed in the 
six fields nearest to the centre of the Coma cluster in each 
of the 2004 and 2007 observing seasons from the archive. 
We stored them in a two-dimensional array whose columns 
held the light curves of individual stars, and whose rows held 
data derived from individual CCD frames. Patterns of corre- 
lated systematic error were identified and removed using the 
SysRem algorit hm of|Tamuz, Mazeh, fc Zucker] ( |2005[ ), as 
implemented by Collier Cameron et al. (20061. At this stage 



the photometric variances are augmented by an additional 
systematic variance that accounts for the actual scatter in 
the measured fiuxes about the final decorrelated solution for 
each frame. The maximum-likelihood algorithm for comput- 
ing this additional variance is described by [Collier Cameron| 
et al. ( 2006 1. The long-term RMS scatter of the decorrelated 



SuperWASP light curve of a non-variable star is 0.004 mag 
at V = 9.5, degrading to 0.01 mag at V = 12.0 as shown in 
Fig. 3 oflCoUier Cameron et al.|(|2006l. 
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4 LIGHT CURVE ANALYSIS 



4.1 Frequency analysis 

We searched for evidence of quasi-sinusoidal light-curve 
modulation in all stars in the field, using the generalised 
Lomb-Scargle periodogram formulation of [Zechmeister fc| 
|Kurster| ( |2009[ ) to fit an inverse variance-weighted floating- 
mean sinusoid to the light curve over a full observing season. 

Each star's light curve comprises a set of observed mag- 
nitudes ra'i £ {m'l, 7712, m^} with associated variances 
G {cyi,o-2, ...,a%}. We form differential magnitudes 



771, — rUi 



where 

The goodness of fit of this constant-magnitude model to the 
light curve is given by the statistic 



2 



We fit a sinusoidal model 



A cos Luti + B sin ujti + C 



with angular frequency cu to the residuals rrii at the times 
ti of observation. We solve for the coefficients A, B and C 



using the algorithms of Zechmeister & Kiirster ( 2009 1 . The 
light curve amplitude at frequency ui is then 



5771 : 



A^ + B^ 



and its associated variance is 

, A^Y&v{A) + B^^s.r[B) 
Var 6771) = ; ^ — -■ 

A^ + B'^ 

We repeat the process at a set of frequencies correspond- 
ing to periods ranging from 1.11 days to 20 days. At each 
frequency u> we compute the statistic 



p{uj) = 



Xo 



x'H 



For each star, the peak value pbcst of the p statistic is 
achieved at the the frequency at which is minimised. 

For a light curve comprising independent measure- 
ments, 



Gumming, Marcy, & Butler ( 1999[ l estimate 



Prob(p > Pbest) = (1 - Pbest 



(iV-3)/2 



(1) 



The false-alarm probability (FAF) for Ai independent fre- 
quency measurements is then 



FAP = l-[Prob(p>pbcst)]' 



(2) 



where M ~ T{f2 — fi) for a data train of total duration 



T searched in the frequency range /i < / < /2 (Gumming 



2004) 



The false-alarm probability is sensitive to the number 
A'^ of independent observations. As |Collier Gameron et al.| 
( 2006 1 found, some forms of correlated error repeat at in- 



tervals of one day, and are not adequately removed by the 
SysRem algorithm. A small misalignment of SuperWASP's 
polar axis causes every stellar image to drift across the frame 
by a few tens of pixels during each night. The vignetting 
pattern of the camera lenses produces discontinuities in the 




Log FAP 



Figure 1. Amplitude of best-fit sinusoidal model as a function of 
false-alarm probability for stars in the field SW1216-I-3126 in the 
2004 season. Stars to the right of the selection threshold set at 
FAP= 10~^ are considered to be candidate rotational variables. 



focal-plane illumination gradient near the edge of the frame, 
which change slightly from night to night and even during 
the course of a single night. The resulting correlated errors 
serve to reduce the number of independent observations from 
A'^ to a lower value A^off- 

We estimated A'ofi by randomising the order of the list 
of dates on which data were obtained. The shuffling proce- 
dure shifts each night's observations in their entirety to a 
new date. This procedure effectively destroys coherent sig- 
nals with periods longer than 1 day, but retains the global 
form of the window function and the effects of correlated 
noise. If many such trials are made, we expect pbost in their 
individual periodograms to exceed the value found for a sin- 
gle randomly-chosen trial roughly half the time. We thus 
estimate a false-alarm probablility FAPo = 0.5 for the the 
strongest peak pbost in the periodogram of the shuffled data. 
Using this estimate, we combine and invert equations [l] and 
[2] to obtain an estimate of the effective number of indepen- 
dent observations in the presence of correlated errors: 



Afci 



31n(l-pbc.t) + 21n(l- (l-fAPo 
ln(l — pbcst) 



.l/M^ 



(3) 



In practice we found that typically A^off — A'^/10. Since Su- 
perWASP observes each field at intervals of 8 minutes or so, 
the effective correlation time is of order 1.5 hours. 

In Fig. [1] we plot the best-fitting modulation ampli- 
tude against false-alarm probability for all stars in the field 
SW1216-)-3126 observed in the 2004 season. We select for 
further study all those light curves showing evidence of pe- 
riodic modulation with false-alarm probabilities less than 
lO"'^. These stars exhibit light curve amplitudes in the range 
0.01 to 0.1 mag, which is typical of rotational variables with 
periods of order a few days. 
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Figure 2. Proper-motion diagram of stars with significant 
(FAP< 10^"^) quasi-sinusoidal variability and periods between 
1.1 and 20 days, in the WASP fields surveyed in 2004 and 2007. 
The red cross denotes the centroid of the field distribution for 
stars with total proper motions less than 25 mas/yr, while the 
red circle marks the mean proper motion of the cluster members. 
The 30 stars with combined kinematic and photometric member- 
ship probabilities greater than 0.5 (see also Table|2j are encircled. 



Figure 3. 2MASS {K, J — K) colour-magnitude diagram of stars 
with significant (FAP< 10~^) quasi-sinusoidal variability and pe- 
riods between 1.1 and 20 days, in the WASP fields surveyed in 
2004 and 2007. Previously-known cluster members detected as ro- 
tational variables and satisfying our cluster membership criteria 
are encircled. The solid curve is the 520 Myr isochrone of |Baraffe| 
|et al.| ( |1998^ for stars of solar metallicity at the distance of the 
Coma Berenices cluster. 



5 CANDIDATE SELECTION 

Having identified a large number (A'^tot = 1613) of rota- 
tional variables in the vicinity of the cluster, the next step 
is to determine which among them are likely to be cluster 
members. We followed the method of Girard et al. (19891, 



defining probability density functions for the proper motions 
of Nc cluster members and Nj = A^;ot — Nc field stars. 

We obtained proper motions and magnitudes in the V 
and 2MASS J and Ks bands for all candidate rotational 
variables by cross-matching with the NOMAD-1 catalogue 



JZacharias et 
iDobbic (20061 



al. 20051. Following Casewell 



we adopt /ic 
9.16 ±0.15 mas yr" 



Jameson, 



& 



= -11.21 ± 0.26 mas yr"^ and 
^ as the mean proper motion of 



the known cluster members, as determined by [van Leeuwen| 
(19991. The proper-motion diagram for the full sample is 



shown in Fig. |2] The field distribution for stars with total 
proper motions less than 25 mas y~^ is approximated by a 
two-dimensional gaussian probability density function hav- 
ing the form 



27rlj f rv5] 



■ exp 



2S2 



9T2 



(4) 



with mean proper motion fj,fa ~ —4.42 mas yr~^, 
—4.39 mas yr~^ and dispersion E/q, — 9.90 mas yr~^ 
8.18 mas yr~^. 

For cluster members the density function depends pri- 
marily on the uncertainties in the proper-motion compo- 
nents aai and asi for the individual stars, since the intrin- 
sic spread in the proper motions of the cluster members is 
smaller than the measurement errors: 




Rank 



Figure 4. Ranked cluster membership probabilities of the top 
50 stars whose light curves recorded by one or more WASP 
cameras m 2004 or 2007 showed significant (FAP< 10-^) quasi- 
sinusoidal variability and periods between 1.1 and 20 days. The 
membership probability is determined from the proper motions 
and K magnitude distributions for cluster and field stars, as de- 
scribed in the text. 30 stars have cluster membership probabilities 
greater than 0.5. 
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2-KOcaOcS 



exp 



2al 



2(7? 



(5) 



Following |Girard et al.| ( |l989[ ), the cluster membership prob- 
abihty for an individual star is 



fc 



fc + ff 



We generalised this approach to take into account the 
magnitude distribution of field stars and proximity to the 
cluster main sequence in a colour-magnitude diagram con- 
structed from the 2MASS K magnitudes and J — K colour 
index (Fig.jsJ. To select likely cluster members using their 
locations in this C-M diagram, we adopted the 520 Myr 
isochrone from Baraffe's (1998) NEXTGEN stellar evolu- 
tion models, for low-mass stars of solar metallicity. The the- 
oretical absolute K magnitudes were converted to apparent 
magnitudes us ing the HIPPARCO S distance of 89.9±2.1 pc 
to the cluster ( van Leeuwen|1999 1 . The Baraffe J ~ K colour 
and K magnitude were converted from the CIT system for 
which they were computed to their 2MASS equivalents us- 
ing the transformation of Carpenter ( 2001| ). We allow an 
intrinsic scatter Scif = 0.15 mag about the theoretical main 
sequence, assuming the depth of the cluster to be of order the 

5 to 6 pc tidal radius estimated by [Qdenkirchen, Soubiran^ 

6 Colin ( 1998 1. The probability density function is assumed 



to be gaussian. 

Because the WASP bandpass is in the optical range, 
the colour-magnitude diagram of the field stars shows Ks 
magnitude to be correlated with J — K colour. The offset in 
K magnitude of a star i from the mean relation is given by 

SKf^ = K,- 11.700 + 2.808((J - K)^ - 0.461), 

with an approximately gaussian intrinsic scatter E/k = 1.38 
mag. 

We multiply the probability densities for the Ks magni- 
tude offset from the field distributions with the field proper- 
motion probability density from Eq. [4] above to obtain 



9S = 



exp 



5K% 

^fK 



V27rE/K 

and similarly for the cluster probability density we obtain 

fc „ „ (^^1 ^^isochronc ) 



(6) 



■ exp 



(7) 



The combined membership probability for an individual star 
is then 



P = 



gc + gf 

In Fig. [5] we plot the cluster membership probabilities 
of the 50 rotational variables that are the most likely cluster 
members, in descending order of rank. The properties of 
the 30 stars with membership probabilities greater than 0.5 
are listed in Table[2] together with cross-identifications with 
previous catalogues of cluster members. 



6 PERIOD-COLOUR RELATION 

The rotation periods of all stars with false-alarm probabili- 
ties less than 0.01 and membership probabilities greater than 
0.5 are plotted as a function of J—K colour in Fig. [5] Among 




Figure 5. Photometric rotation period as a function of 2MASS 
(J — K) colour for stars with significant (FAP< 10"^) quasi- 
sinusoidal variability and cluster membership probability greater 
than 0.5, in the WASP fields surveyed in 2004 and 2007. Periods 
measured for candidates with cluster membership probabilities 
greater than 0.85 are denoted by filled circles; those with mem- 
bership probabilities between 0.5 and 0.85 are denoted by open 
circl es. Additional stars from [Radick, SkiflF, &: Lockwood| | |1990[ l 
and |Marilli, Catalano, fc Frascaj l |1997| are denoted by '-|-' and 
'x' symbols respectively. 



these, the periods of stars with cluster membership proba- 
bilities greater than 0.85 are denoted by solid symbols. They 
follow a tight period-colour relation. The properties of these 
stars are listed in Table [2] including cross-identifications 
with previous catalogues of cluster members. The rotational 
parameters of individual light curves that yielded significant 
detections in the 2004 and 2007 seasons are listed in Table El 

Several light curves yield periods that follow a sec- 
ond sequence in Fig. [5] with periods close to half those 
on the main relation. A similar secondary period-colour se- 
quence was noted by [Hartman et al.| ( |2008[ l in their re- 
cent study of the comparably-aged cluster M37. Several of 
them were recorded independently in two or more of the 
WASP cameras, and some in different seasons, and were 
found to jump between the two relations. The K4 dwarf 
ISWASP J121135. 15+292244.5, for example, yields P = 
10.37 d in 2004, and 5.07d in 2007. Similarly, the Kl dwarf 
J122856.43+263257.4 yielded periods of 4.79 d in 2004 and 
9.26d in 2007. The K2 dwarf ISWASP J121253.23-f 261501.2 
yields 10.38d in 2004 and 5.44d in 2007, though only the lat- 
ter measurement has FAP< 0.01 and is listed here. 

These stars almost certainly rotate with periods near 
the main relation. The reason for the intermittent frequency 
doubling is apparent from Fig. [6] The 2004 light curve yields 
a period half that of the 2007 light curve, but the longer pe- 
riod is clearly present in both seasons. Evidently the appar- 
ent doubling of these stars' rotational frequencies at some 
epochs of observation arises from two dominant starspot 
groups on opposing stellar hemispheres giving a double- 
humped light curve. We conclude that the true periods of the 
light curves on the lower sequence are twice those derived 
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1SWASP Jl 22856.43+2632 P= 4.79d iV = 0.006 



1SWASP J 122856.43+2632 FAP=0.000094 i)(^= -1076.6 
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Figure 6. Light-curves and pcriodograms for ISWASP J122856. 43+263257. 4 observed in 2004 (top) and 2007 (bottom). Eacii set of four 
panels stiows the unphased liglit-curve at upper left. Occasional gaps in the data were caused by inclement weather. The periodogram 
of p = Ax^/Xo versus frequency is at upper right, with horizontal dashed lines denoting thresholds of false-alarm probability 0.1, 0.01 
and 0.001. The vertical bar indicates the true rotational frequency. The light curve phased on the best-fit period for the year concerned 
is shown at lower left, and the sine-wave amplitude periodogram at lower right. The 2004 data illustrate the frequency-doubling that 
occurs when major spot complexes are present on opposite stellar hemispheres, giving a double-humped light curve. 
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Table 2. Cross-identifications, 2MASS photometry and proper motions of candidate rotational variables with membership probabilities 
greater than 0.5. The ISWASP identifier encodes the J2000.0 coordinates of the object, d denotes the angular distance of the star from 
the cluster centre. The spectral types are estimated from 2MASS photometry using the calibration of |Blackwell fc Lynas-Gray| ( |1994| . 
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from sine fitting. After correcting for frequency-doubling, 
the period-colour relation is as shown in Fig. [7] 

Among the stars with membership probabilities greater 
than 0.85, only one outlier is seen with a period that departs 
significantly from the main relation. This star, ISWASP 
J124930.43+253211.1 {J - K = 0.682), exhibited a clear 
period of 19.51d in 2004, with a weak secondary minimum. 
In 2007 the period appears to be 9.41 d, but the minima 
alternate in depth indicating the true period to be 18.82d. 



Meibom, Mathieu, & Stassun ( 2009 1 noted the presence of 



several similarly anomalous slow rotators in their rotation 
study of M35, and proposed that partial or complete tidal 
synchronisation in stellar binary systems could be responsi- 
ble. Radial- velocity observations of this star would be desir- 
able to check whether it really belongs to the cluster, and 
to test for evidence of binarity. 

Half of the 23 rotational variables with membership 
probabilities greater than 0.85 lie within 4.2 degrees of the 
cluster centre, while the other half lie in the periphery of 
the cluster, between 4 and 11.5 degrees from the cluster cen- 
tre. Among them, 14 appear to be previously uncatalogued 
cluster members. Rotational variability thus appears to be 



a useful method for identifying new members of the cluster 
and its extended halo. 

In Figs.[5]and[7]we also show previously-published pho- 
tometric periods for Melotte 111 stars T65, T76, T85 and 



T132 ( Radick, Skiff, fc Lockwood|1990|) and T203, T213 and 
T221 ( Marilh, Catalano, fc Frasca|1997[ ), with their 2MASS 
colours. 

The photometric periods found by |Marilli, Catalario^ 



fc Frascal ( [T997| ) for T203 and T221 lie close to the Super- 
WASP period-colour relation. The 4.7-day period they found 
for T213 is very close to half the period predicted by Eq. |8] 
The sine-fitting method reveals a 9.51-day modulation in the 
WASP hght curve of this star, but the false-alarm probabil- 
ity was below the threshold value 0.01 that would qualify it 
as a secure detection. We adopt a period of 9.4 days for this 
object, twice that reported by Marilli et al. 

Together with the four Radick stars and the three Mar- 
illi stars, but excluding the anomalous slow rotator ISWASP 
J124930.43+253211.1, the WASP hght curves for stars with 
high cluster-membership probabilities and low false-alarm 
probabilities lie along a well-defined, almost linear period- 
colour relation. A linear least-squares fit to the period-colour 
relation of these stars yields 
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Table 3. Sine-fitting periods for light curves of candidate rotational variables in the vicinity of Melotte 111, with cluster membership 
probabilities greater than 0.5. The year of observation is listed, with a designator incorporating the RA and Dec of the field centre and 
a three-digit camera identifier for each light curve. The number N^y^^ of observations is listed, followed by the period of the peak in the 
periodogram yielding the greatest improvement in the statistic. Columns 6 and 7 give the amplitude of the fitted sinusoid and the 
false-alarm probability of the signal detection. Periods suspected of being half the true period appear in parentheses. 
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43-1-263257 


4 


2004 


SW1244+2427 


.101 


2268 


(4.79) 





006 





001287 


J122856 


43-1-263257 


4 


2007 


SW1242-I-2418 


.141 


3833 


9.26 





007 





006194 


J122942 


15-1-283714 


6 


2007 


SW1222-I-3000 


.144 


2894 


16.05 





004 





003018 


J123231 


07-1-351952 


2 


2004 


SW1243+3126 


.102 


2259 
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013 





000240 


J123941 


99-1-213458 





2004 


SW1244+2427 


.101 


2264 
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P = 9.71 + 10.68( J - 0.528). 



(8) 



The RMS scatter of the observed periods about the fit- 
ted period-colour relation is 0.19 days - just 2 percent of 
the mean period - and is not improved by the inclusion of 
quadratic or higher terms. This scatter is so small that much 
of it is attributable to year-to-year changes in the photomet- 
ric modulation periods of the cluster stars themselves. There 
are seven instances in which periods were measured for the 



same star in both 2004 and 2007. These season-to-season 
period differences are 0.09, 0.11, 0.16, 0.23, 0.32 and 0.49 
days. The season-to-season differences in the period deter- 
minations for a given star average 0.23 days, comparable to 
the scatter about the mean relation. Radick et al. ( 1995 1 



reported finding very similar season-to-season period varia- 
tions of between 2 and 8 percent in their 12-year study of 
rotation in the Hyades. 
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Figure 7. As for Fig. [s] with the periods of the light curves on 
the lower branch having been doubled. The solid line is a least- 
squares linear fit to the period-colour relation, and the dashed 
lines delineate the scatter about the mean relation expected from 
cyclic changes in spot latitude on differentially-rotating stars. 
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Figure 8. Period-colour relations for Coma Berenices (filled cir- 
cles) and the Hyades (open circles). The Hyades J — K colours 
are taken from 2MASS, and the Hyades rotation periods from 
IRadick et al] | |1987| |. 



There are two possible explanations for this trend: lat- 
itudinal diflerential rotation and core-envelope decoupling, 
both of whose effects will be superimposed on a spread at 
all spectral types resulting from a range of disc-locking life- 
times. 

Differential rotation on the stars themselves will lead 
to secular changes in the modulation period, if activity cy- 
cles are present which cause the main spot belts to drift in 
latitude. Doppler imaging and line-profile studies of differ- 
ential rotation in young main-sequence stars show a strong 
decrease in latitudinal differential rotation along the main 
sequence. F and G dwarfs exhibit stronger shear than mid- 



to-late K dwarfs which rotate almost as solid bodies ( Barnes 
et al.|2005[|Reiners|200"6l ). |Collier Cameron| ( |2007[ ) gives the 
empirical relation for the equator-to-pole difference in rota- 
tional frequency 



An = 0.053(rgfj/5130)* ^ad d"\ 



(9) 



We expect the resulting scatter in rotation frequency 
to be no more than half this amount, if the main active 
belts migrate across a similar range of latitudes as the solar 
butterfly diagram occupies on the Sun. The corresponding 
range of rotation periods at each colour is shown in Fig. [7] 
The dispersions about the main period-colour reported here 



and by Hartman et al. ( 2008 1 in M37 are comparable to the 
scatter expected from differential rotation for the F and G 
stars. The K stars also show season-to-season period changes 
that are most easily explained by a degree of differential ro- 
tation somewhat greater than Doppler imaging studies have 
found among more rapidly-rotating stars of similar mass. 
At all colours we find the season-to-season period changes 
to contribute a significant fraction, though perhaps not all, 
of the overall scatter about the mean period-colour relation. 

Additional residual scatter is expected at all spectral 
types if the early spin rates of stars are regulated by torque 
balance with accretion discs having a range of lifetimes. Such 



"disc-locking" models are reasonably successful at explain- 
ing the distribution of stellar spin rates on the zero-age main 
sequence ([Collier Cameron, Campbell, fc Quaintrell 1995 



Keppens, MacGregor, fc Charbonneau|1995||Krishnamurthi| 
et al.|1997 ,, AUain 1998). 



Once the star decouples from the disc, angular momen- 
tum is lost through the stellar wind alone, but the spindown 
behaviour of the stellar surface depends critically on the de- 
gree of coupling between the star's radiative interior and 



convective envelope. Stauffer & Hartmann (19871 sought to 



explain the rapid spindown of G dwarfs between the ages of 
the 50 Myr-old a Per cluster and the 70 Myr-old Pleiades by 
suggesting that magnetic braking early in a star's life may 
spin down only the outer convective layers at first, leaving 
a decoupled reservoir of angular momentum in the stellar 
radiative core. If the coupling timescale on which the stored 
angular momentum is transported from the core back into 
the envelope is longer than a few tens of Myr, outward an- 
gular momentum transport will spin up the envelope at in- 
termediate ages. This will counter the effects of magnetic 
braking and delay full convergence to a single period-colour 
relation. The effect would be lessened in K dwarfs, whose 
radiative interiors comprise a smaller fraction of the total 
stellar moment of inertia. 

Several authors have modelled and discussed the impli- 
cations of core-envelope decoupling for cluster period distri- 
butions at various ages (iLi fc ColUer Cameron|1993l |Collier| 
[Cameron fc Lill994 Krishnar nurthi et al.|1997||Allain|l'998| 
"Wolff, Strom, & Hillenbrand 2004'; 'Soderbl om, Jones, fc Fis-| 



|cher,,2001^ . The models in these papers indicate that the 
weakest core-envelope coupling that can ensure the observed 
near solid-body rotation in the solar interior at an age of 4.6 
Gyr would produce a dispersion by a factor two or more in 
the spin rates of solar-type stars at the age of the Hyades 
and Coma. The observed scatter among Coma stars with 
mid-F to mid-G spectral types is, however, only of order 2 
percent. Since much of this is demonstrably caused by differ- 
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ential rotation, the evolutionary contribution to the scatter 
must be small. This argues strongly for coupling timescales 
substantially shorter than the cluster age. 



7 GYROCHRONOLOGICAL AGES OF COMA 
AND HYADES 

The tight period-colour relation for the rotational variables 
found within 12 degrees of the core of the Coma Berenices 
cluster, sharing the proper motion and lying close to the 
cluster main sequence, contrasts strongly with the bimodal 
distribution of rotation periods seen in younger clusters. The 
absence of a population of rapid rotators in our survey is 
not, however, surprising because our faint limit extends only 
to the boundary between K and M spectral types. In the 
Hyades, rapid rotation is only seen at spectral types later 
than our faint limit. At earlier spectral types, the period- 
colour relation in Coma resembles the pattern of rotation 
rates observed in the Hya des by|Radick et al.] ( |1987[ ), and 



subsequently updated by Radick et al. (19951. In Fig [s] 



we show the period-colour relations for both clusters, us- 
ing 2MASS J — K colours for the Hyades stars. All but five 
of the Radick et al. stars were included in the Hipparcos 



survey of the Hyades by Ferryman et al. (19981, and were 
identified as cluster members from their parallaxes, proper 
motions and radial velocities. 

The relative rotational ages of the two clusters can 
be inferred from their respective period-colour relations us- 
ing gyrochronology, assuming a simple rotational spindown 
model. We follow a methodology similar to that of IBarnes] 
( |2007[ ), assuming the rotation period of a given star to be 
a separable function of age and colour. We divide the pe- 
riods of stars in both clusters by the periods derived from 
their J ~ K colours at the age of Coma using eq. |8] Assum- 
ing that rotation period increases as the inverse square root 
of age, the age of an individual star relative to the fiducial 
rotational age tcoma of the Coma population is given by 



t = tCon 



9.71 -I- 10.68(J - A" - 0.528) 



(10) 



The relative gyrochronological ages of individual stars 
in the Hyades and Coma are shown as a function of J — K in 
Fig. |9] The sample means for the relative ages of the Coma 
and Hyades stars are 1.002±0.015 and 1.058±0.065 respec- 
tively. These values exclude the lone Hyades rapid rotator 
with J — K — 0.77, and the anomalous slow rotator dis- 
cussed above at J ~ K = 0.682. The difference between the 
gyrochronological ages of the two clusters is not statistically 
significant, owing to the high RMS scatter in the relative 
age determinations for the Hyades stars. The distribution 
of ages for individual Coma stars is much tighter, showing 
an RMS scatter of only 0.090 in relative age. If we adopt 



625 Myr as the age of the Hyades cluster (Ferryman et al 



|1998[ ), the relative gyrochronological age of Coma is found 
to be 590.7±40.9 Myr. The uncertainty is dominated by the 
scatter of the Hyades measurements, as seen from the cumu- 
lative distributions in Fig. |10| The mean gyrochronological 
age of the Coma sample is thus 34 ± 41 Myr younger than 
that of the Hyades sample. This is in reasonable agreement 
with the recent comparative study of the isochrone ages of 
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Figure 9. Gyrochronological ages for individual light curves of 
stars in Coma Berenices (filled circles) and the Hyades (open 
circles) , derived from Eq. |10| 
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Figure 10. Cumulative distributions of gyrochronological ages 
for individual light curves of stars in Coma Berenices (filled cir- 
cles) and the Hyades (open circles), derived from Eq. |10| The 
mean age of the Coma Berenices stars is 591 Myr; that of the 
Hyades is 625 Myr. 



to be up to 100 Myr younger than the Hyades. Adopting 
^Coma = 591 Myr as the age of Coma, we derive the period- 
age-colour relation. 



P = (9.71 + 10.68( J -K- 0.528)) ^7^591 Myr. 



(11) 



Eq. [TT] yields a predicted sidereal solar photometric 
modulation period of 21.8 days for an assumed (J — isT)© = 
0.354 ( [Ramirez fc Melendez|2005| l and a solar age = 4.56 
Gyr. This is shorter than the observed 26.1-day sidereal pe- 



these clusters by King & Schuler ( 2005 1 , who found Coma 



riod of sunspots in the mid-latitude active belt (Donahue, 
|Saar, fc Baliunas]|1996[ ). Conversely, the Sun's 26.1-day pe- 
riod yields an age of 6.5 Gyr using Eq. [lOlwith fcoma = 591 
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Table 4. J-K colours and spin periods for the Sun and old main-sequence stars with measured rotation periods and asteroseismological 
age determinations. The J — K colours were determined using the effective temperatures and interpolating formulae of Ramirez ^ 
Melendez] ||2005 I. The rotation periods are from [Donahue, Saar, &: Baliunas] ( |1996^ for the Sun, |Barnes| l |2007| for a Cen, and from 
Stimets fc Gilesll ligSO] and |Noyes et al.| ( |1984] for 70 Oph A. The last row is based on the supposition that the true period of 70 Oph A 



is twice the measured period. The asteroseismological ages for a Cen and 70 Oph A arc from Eggenberger et al. (2004, 2008). The gyro 
ages are estimated for two different values of the magnetic braking power-law index i 



Star J — K Period Gyro Age Gyro Age Seismo 

(days) (6 = 0.5) {b = 0.56) Age (Gyr) 



Sun 


0.354 


26.1 


± 


2.0 


6.5 


± 


1.0 


5.0 ± 


0.7 


4.57 


a Cen A 


0.35 


28.0 




3.0 


7.6 


± 


1.7 


5.8 ± 


1.1 


6.5 ±0.3 


a Cen B 


0.49 


36.9 




1.8 


9.3 


± 


0.9 


6.9 ± 


0.6 


6.5 ±0.3 


70 Oph A 


0.55 


19.9 




0.5 


2.4 


± 


0.1 


2.0 ± 


0.1 


6.2 ± 1.0 


70 Oph A 


0.55 


39.8 


± 


1.0 


9.5 


± 


0.5 


7.0 ± 


0.3 


6.2 ± 1.0 



Myr. While an idealised braking law of this kind is conve- 
nient, it is unlikely to be appropriate for real stars whose 
moments of inertia evolve throughout their main-sequence 
lifetimes. A convenient first-order correction may be a simple 
adjustment to the index b of the braking power law P oc , 
as Barnes (20071 has proposed. 

The major obstacle to fine-tuning the braking power- 
law index is the lack of old stars with reliably-determined 
ages and spin periods. A very few stars have already had 
their spin periods determined from the rotational modula- 
tion of their chromospheric Ca II H&iK emission, and their 
ages determined from their solar-like p-mode oscillation pat- 
terns, namely a Cen A and B, and 70 Oph A. [Barnes] ( |2007[ ) 
found a Cen A and B to have comparable gyrochronologi- 
cal ages, from their respective spin periods of 28.0 and 36.9 
days; [Eggenberger et ah] ( [2004| ) derived an asteroseismolog- 
ical age of 6.5 ± 0.3 Gyr for the a Cen system. The spin 
period of 70 Oph A was found by [Stimets fc Giles| ( [1980[ ) to 
be 20.1 days and by [Noyes et al.[ ( [1984[ ) to be 19.7 days from 
the same set of Mt Wilson Ca II H&cK modulation data. 



Eggenberger et al.| ( 2008 1 determine an asteroseismological 
age of 6.2 ± 1.0 Gyr for 70 Oph A. 

In Table [4[ we show that the rotational ages of these 
stars are best reconciled with the asteroseismological ages 
by adjusting the braking power law index from b — 0.5 to 
b — 0.56. While neither value of the index b gives an age 
for 70 Oph A that agrees with the asteroseismological age 
of 6.2 Gyr at a spin period of 20 days, we note that Ca II 
HSzK light curves can also exhibit double-peaked modula- 
tion in some seasons. It is possible that the true spin period 
of 70 Oph A is close to 40 days; future analyses of rota- 
tional splitting in the low-order p-modes of 70 Oph A would 
provide a valuable independent check on this star's rotation 
rate. If the longer period is correct, the asteroseismological 
and gyro ages can be reconciled using the same 6 = 0.56 that 
fits the Sun and both components of a Cen. The best-fitting 
period-colour-age relation is then 



t = 591 



9.30 



P 

10.39(J- 



K - 0.504) 



1/0.56 



Myr. (12) 



A departure of the braking power-law index from the 
b = 0.5 of the Skumanich relation also seems to be needed 
in the comparison between K dwarfs in M35 and the Hyades 



( Meibom, Mathieu, fc Stassun|20"09| ). [IIartman et al.| ( [2008 l 
also find that they cannot fit M37, Hyades and the Sun 
simultaneously with a braking law that yields an asyrnptotic 



-1/2 



age dependence for the stellar spin rate. 



Barnes 



(20071 



and Mamajek fc Hillenbrand] (2008) find similar departures 
from the Skumanich relation in their gyrochronological age 
calibrations for cluster and field stars. For the purposes of 
the present investigation, however, the ages of the Hyades 
and Coma are in any case so similar that such small changes 
in the braking index will have little effect on the ratio of the 
two clusters' ages. 



8 CONCLUSIONS 

We have shown that inverse variance-weighted sine-fitting 
searches yield reliable period determinations from Super- 
WASP data in stars as faint as Ks = 10, with amplitudes of 
rotational variability as low as 0.003 mag and periods of or- 
der 10 days recovered reliably with generalised Lomb-Scargle 
false-alarm probabilities less than 0.01. This has enabled us 
to carry out the first comprehensive study of rotation among 
the F, G and K dwarfs in the Coma Berenices open cluster. 
We find a narrow, nearly linear relation between rotation 
period and colour from mid-F to late-K spectral types. 

The scatter of the individual stellar rotation rates about 
the mean period-colour relation is surprisingly small, as has 



previously been noted by Radick et al. (19871 from their 
study of the Hyades, by [Hartman et aTf j 2008P and [Messina[ 



et al. ( 2008 1 from their more recent studies of M37, and 



by Meibom, Mathieu, fc Stassun (20091 from their study 



of M35. Indeed, our period-colour relation for Coma has a 
markedly lower scatter about the mean relation than Radick 
et al. found for the Hyades. It is not clear whether the dif- 
ference in scatter is intrinsic, or simply refiects the much 
denser sampling of the SuperWASP observations. Some of 
the rotation periods in the Radick et al. Hyades sample were 
derived from as few as 30 observations per star obtained over 
an interval of 5 months. The SuperWASP periods in Coma 
are derived from between 1000 and 3800 observations over 
a four-month period in each observing season. 

The RMS scatter about the mean period-colour rela- 
tion in Coma is less than 5 percent. A substantial fraction 
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of this scatter appears to be attributable to season-to-season 
changes in period, caused by differential rotation and secu- 
lar changes in active-region latitudes. The remaining few- 
percent residuals about the fitted period-colour relation are 
consistent with the not-quite-complete convergence in spin 
rates expected from rotational models at this age. We do 
not see the large scatter expected from models with core- 
envelope coupling timescales longer than about fOO Myr. 
Even solid-body spindown models leave a residual scatter 
in rotation rates at the age of Coma and the Hyades. This 
residual scatter depends somewhat on the detailed form of 
the braking law among the fastest rotators at early ages, 
but is seldom less than about 10 percent in the models cited 
above. 

Although the stars in these clusters have not entirely 
forgotten their original spin rates, the degree of convergence 
is good enough to yield reasonably accurate gyrochronolog- 
ical age estimates. The scatter in the period-colour relation 
propagates into the distribution of ages for individual cluster 
stars derived from a linear pe riod-colour relatio n assuming 
a simple t~^^'^ spindown law (Skumanich 19721. The most 
remarkable feature of Fig. [10 is the sharpness of the age 
distribution inferred for Coma. The RMS scatter in the gy- 
rochronological ages for the Coma stars is only just over 9 
percent of the mean age derived for the ensemble. This high 
degree of internal consistency is echoed in M37, M35 and 
the Hyades, as other authors have noted previously. 

We have used J — if as the colour index of choice for 
defining a period-colour-age relation based on the Coma 
data, partly because of the ready availability of well- 
determined JHK colours in the 2MASS point-source cat- 
alogue for stars in the SuperWASP magnitude range, and 
partly because J — K is less subject to reddening and metal- 
licity effects than the more widely-used B — V (Houdashelt, 
Bell, fc Sweigart|2000 1. Differences in metallicity from clus- 
ter to cluster may affect the B — V colour at a given mass, 
and the period to which stars of a given mass converge at a 
given age. We raise the latter possibility because the angular 
momentum loss rate in a thermally- driven wind depends on 
the wind temperature as well as the total open magnetic flux 
threading the wind ( Mestel fc Spruit|1987 Collier Cameron 



fc Li 19941. A dependence of wind temperature on metal- 



licity would lead to systematic differences in the asymptotic 
rotation rate for stars of a given mass and age in clusters 
with different abundances. 

Despite these caveats, the relative ages derived here for 
Coma and the Hyades are in close agreement with the dif- 
ference in their isochrone ages, even though, according to 



Cayrel de Strobel et al. (19971, the metallicity of Coma is 



sub-solar by about 0.05 dex, while the Hyades is overabun- 
dant by about -1-0.15 dex. Our results thus support the con- 



tention of Barnes (20071 that a properly-calibrated stellar 



rotational clock can provide a viable means of age determi- 
nation in F, G and K stars that rotate more slowly than 
their counterparts of the same colour in the Hyades, Coma 
Berenices and M37 clusters. 

Some fine-tuning will still be needed to obtain gy- 
rochronological ages of the best possible precision for main- 
sequence F, G and K stars. We contend that asteroseis- 
mology offers a promising way forward, providing precise 
measurements of both stellar ages and spin periods in ma- 
ture main-sequence stars. Although such analysis is exacting 



and only feasible for very bright stars, it will provide valu- 
able gyrochronological calibration standards at later ages. 
Gyrochronological ages can then be determined for much 
larger numbers of fainter stars, in a manner analogous to 
the cosmological distance ladder. 
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